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^ ■ ABSTRACT 

^ ■ Interactions of heavy Majorana neutrinos in the thermal phase of the early uni- 

verse may be the origin of the cosmological matter-antimatter asymmetry. Suc- 
V) . cessful baryogenesis, independent of initial conditions, is possible for neutrino 

masses in the range 10 -3 eV < m, < 0.1 eV. Remarkably, this mass window is 
r ^ ■ consistent with the evidence for neutrino masses from oscillations. 

t> 
o , 

^ ■ 1. Matter- Antimatter Asymmetry 

m 

One of the main successes of the standard early-universe cosmology is the pre- 
diction of the abundances of the light elements, D, 3 He, 4 He and 7 Li. Agreement 
between theory and observation is obtained for a certain range of the parameter t)b, 
the ratio of baryon density and photon densitj^ , 
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V BBN = ^B = ( 2 q _ Q 2) x 10 -10 ? (1) 



n 7 



x. 

where the present number density of photons is n 7 ~ 400/cm 3 . Since no significant 
amount of antimatter is observed in the universe, the baryon density coincides with 
the cosmological baryon asymmetry, r] B = {n B — n§)/n 7 . 

The precision of measurements of the baryon asymmetry has dramatically im- 
proved with the observation of the acoustic peaks in the cosmic microwave back- 
ground radiation (CMB). Most recently, the WMAP Collaboration has measured the 
baryon asymmetry with a (Iff) standard error of ~ 

V°b MB = (6.1183) x 10- 10 . (2) 

A matter-antimatter asymmetry can be dynamically generated in an expanding 
universe if the particle interactions and the cosmological evolution satisfy Sakharov's 
conditions^, 
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• baryon number violation , 

• C and CP violation , 

• deviation from thermal equilibrium . 

Although the baryon asymmetry is just a single number, it provides an important 
relationship between the standard model of cosmology, i.e. the expanding universe 
with Robertson- Walker metric, and the standard model of particle physics as well as 
its extensions. 

At present there exist a number of viable scenarios for baryogenesi They 
can be classified according to the different ways in which Sakharov's conditions are 
realized. In grand unified theories baryon number (B) and lepton number (L) are 
broken by the interactions of gauge bosons and leptoquarks. This is the basis of 
classical GUT baryogenesiJ^ . Analogously, the lepton number violating decays of 
heavy Majorana neutrinos lead to leptogenesi P. In the simplest version of leptoge- 
nesis the initial abundance of the heavy neutrinos is generated by thermal processes. 
Alternatively, heavy neutrinos may be produced in inflaton decays, in the reheating 
process after inflation, or by brane collisions. The observed magnitude of the baryon 
asymmetry can be obtained for realistic neutrino masses. 

The crucial deviation from thermal equilibrium can also be realized in several 
ways. One possibility is a sufficiently strong first-order electroweak phase transition 
which would make electroweak baryogenesis possible. For the classical GUT baryo- 
genesis and for thermal leptogenesis the departure from thermal equilibrium is due 
to the deviation of the number density of the decaying heavy particles from the equi- 
librium number density. How strong this departure from equilibrium is depends on 
the lifetime of the decaying heavy particles and the cosmological evolution. 

A crucial ingredient of baryogenesis is the connection between baryon number and 
lepton number in the high-temperature, symmetric phase of the standard model. Due 
to the chiral nature of the weak interactions B and L are not conserve At zero 
temperature this has no observable effect due to the smallness of the weak coupling. 
However, as the temperature approaches the critical temperature T c of the electroweak 
phase transition, B and L violating processes come into thermal equilibrium^. 

The rate of these processes is related to the free energy of sphaleron-type field 
configurations which carry topological charge. In the standard model they lead to an 
effective interaction of all left-handed fermionJ^ (cf. fig. Q), 




(3) 



which violates baryon and lepton number by three units, 



AB = AL = 3 . 



(4) 
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Figure 1: One of the 12-fermion processes which are in thermal equilibrium in the high-temperature 
phase of the standard model. 



The sphaleron transition rate in the symmetric high-temperature phase has been eval- 
uated by combining an analytical resummation with numerical lattice techniques^. 
The result is, in accord with previous estimates, that B and L violating processes are 
in thermal equilibrium for temperatures in the range 

T EW ~ 100 GeV <T < T SPH ~ 10 12 GeV . (5) 

Sphaleron processes have a profound effect on the generation of the cosmological 
baryon asymmetry. Eq. 0] suggests that any B + L asymmetry generated before the 
electroweak phase transition, i.e., at temperatures T > T EW , will be washed out. 
However, since only left-handed fields couple to sphalerons, a non-zero value of B + L 
can persist in the high-temperature, symmetric phase if there exists a non-vanishing 
B — L asymmetry. An analysis of the chemical potentials of all particle species in the 
high-temperature phase yields the following relation between the baryon asymmetry 
and the corresponding L and B — L asymmetries, 

(B) T = c s (B - L)t = -^-r(L) T ■ (6) 

c s - 1 

Here cs is a number 0(1). In the standard model with three generations and one 
Higgs doublet one has c s = 28/79. 

An important ingredient in the theory of baryogenesis is also the nature of the 
electroweak transition. A first-order phase transition yields a departure from ther- 
mal equilibrium. Since in the standard model baryon number, C and CP are not 
conserved, it is conceivable that the cosmological baryon asymmetry has been gen- 
erated at the electroweak phase transition. Detailed studies during the past years 
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Figure 2: Effective, lepton number violating interaction. 



have shown that for Higgs masses above the present LEP bound of 114 GeV elec- 
troweak baryogenesis is not viable, except for some supersymmetric extensions of the 
standard modeff^. In particular, the electroweak transition may have been just a 
smooth crossover, without any departure from thermal equilibrium. In this case it's 
sole effect in the cosmological evolution has been to switch off the B — L changing 
sphaleron processes adiabatically. 

Based on the relation (JHJ) between baryon and lepton number we then conclude 
that B — L violation is needed to explain the cosmological baryon asymmetry if 
baryogenesis took place before the electroweak transition, i.e. at temperatures T > 
Tew ~ 100 GeV. In the standard model, as well as its supersymmetric version and 
its unified extensions based on the gauge group 577(5), 77 — L is a conserved quantity. 
Hence, no baryon asymmetry can be generated dynamically in these models and one 
has to consider extensions with lepton number violation 13 . 

The remnant of lepton number violation at low energies is the appearance of an 
effective AL = 2 interaction between lepton and Higgs fields (cf. fig. EJ), 

£al=2 = \ fij c lidP + n - c - • ( 7 ) 

Such an interaction arises in particular from the exchange of heavy Majorana neu- 
trinos. In the Higgs phase of the standard model, where the Higgs field acquires a 
vacuum expectation value, it gives rise to Majorana masses of the light neutrinos u e , 
Ufj, and v T . 

Lepton number violation appears to be necessary to understand the cosmological 
baryon asymmetry. However, it can only be weak, since otherwise any baryon asym- 
metry would be washed out. The interplay of these conflicting conditions leads to 
important contraints on neutrino properties and on extensions of the standard model 
in general. 

b In the case of Dirac neutrino masses, where the Yukawa couplings of right-handed neutrinos 
are very small, one can construct models where an asymmetry of lepton doublets is accompanied 
by an asym metry o f right-handed neutrinos such that the total lepton number is conserved and 
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2. Leptogenesis 



Lepton number is naturally violated in grand unified theories (GUTs). The uni- 
fication of gauge couplings at high energies suggests that the standard model gauge 
group is part of a larger simple group, 

G SM = £7(1) x SU{2) x SU(3) c SU(5) c 50(10) .... (8) 

The simplest GUT is based on the gauge group SU(5}P. Here quarks and leptons 
are grouped into the multiplets, 

10 = (q L , u c R , e%) , 5* = (d R , l L ) , (1 = u R ) . (9) 

Unlike gauge fields, quarks and leptons are not unified in a single multiplet. In 
particular, right-handed neutrinos are not needed in SU(5) models. Since they are 
singlets with respect to SU(5), they can have Majorana masses M which are not 
controlled by the Higgs mechanism. 

The three SU(5) multiplets can have Yukawa interactions with two Higgs fields, 
H x {b) and # 2 (5*), 

C = huijlOilOjH^) + ^5*40^2(5*) + ta i 5* i l i ffi(5) + M^hlj ■ (10) 

Electroweak symmetry breaking then leads to quark and charged lepton mass matrices 
and to the Dirac neutrino mass matrix mo = h u V\, where v\ = (Hi). The Majorana 
mass term, which violates lepton number (AL = 2), is invariant under SU(5), and the 
Majorana masses can therefore be much larger than the electroweak scale, M 3> v. 

All quarks and lep tons of one generation are unified in a single multiplet in the 
GUT group SO(10^3, 

16 = 10 + 5* + 1. (11) 

Right-handed neutrinos are now required by the fundamental gauge symmetry, and 
the theory contains all ingredients needed to account for the recent evidence for 
neutrino masses and mixings. The seesaw mechanism!^ explains the smallness of 
the light neutrino masses by the largeness of the heavy Majorana masses M. The 
theory predicts six Majorana neutrinos as mass eigenstates, three heavy (jV) and 
three light (z/), 

iV ~ u R + v% : m N ~ M ; (12) 
v ~ v L + v c L : m v = -m D l- m T D . (13) 
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In the simplest pattern of symmetry breaking, B — L, a subgroup of SO(10), is 
broken at the unification scale Lgut- If Yukawa couplings of the third generation are 
(9(1), as it is the case for the top-quark, one finds for the corresponding heavy and 
light neutrino masses: M3 ~ Lgut ~ 10 15 GeV and 7713 ~ v 2 /M% ~ 0.01 eV. It is very 
remarkable that the light neutrino mass m 3 is of the same order as the mass differences 
(Am^) 1 '' 2 and (Am^ tm ) 1//2 inferred from neutrino oscillations. This suggests that, via 
the seesaw mechanism, neutrino masses indeed probe the grand unifiation scale! The 
difference of the observed mixing patterns of quarks and leptons is a puzzle whose 
solution has to be provided by the correct GUT model. Like for quarks and charged 
leptons one expects a mass hierarchy also for the right-handed neutrinos. For instance, 
if their masses scale like the up-quark masses one has Mi ~ 10~ 5 M 3 ~ 10 10 GeV. 

The lightest of the heavy Majorana neutrinos, N\, is ideally suited to generate the 
cosmo logical baryon asymmetry^ 1 . Since it has no standard model gauge interactions 
it can naturally satisfy the out-of-equilibrium condition. Ni decays to lepton-Higgs 
pairs then yield a lepton asymmetry (L)t 7^ 0, which is partially converted to a 
baryon as ymmetry (B)t 7^ 0. The generated asymmetry is proportional to the CP 
asymmetr y^^^ ^ in ^-decays. In the case of the standard model with one Higgs 
doublet, i.e. Hi = H% = <f), the CP asymmetry is conveniently written in the following 
form, 



£1 



T(Ni -> Z0) - T(Ni -> g) 
T{Ni -> l(j>) + T{Ni — U) 
_3 1 i m ( h \ mh *\ n 4 ) 

where the seesaw mass relation (|T3|) has been used. 

From the expression (fTlj) one easily obtains a rough estimate for Ei in terms of 
neutrino masses. Assuming dominance of the largest eigenvalue of m u , phases 0(1) 
and approximate cancellation of Yukawa couplings in numerator and denominator 
one finds, 

3 Mim 3 Mi 

16/T v z M 3 

where we have again used the seesaw relation. Hence, the order of magnitude of the 
CP asymmetry is approximately given by the mass hierarchy of the heavy Majorana 
neutrinos. For M1/M3 ~ 10~ 5 ond has Ei ~ 10~ 6 . 

Given the CP asymmetry S\ one obtains for the baryon asymmetry, 

n B -n B k 9 

t)b = = -zc s ei ~ 10 . (16) 

n 7 / 

Here / ~ 10 2 is the dilution factor which accounts for the increase of the number 
of photons in a comoving volume element between baryogenesis and today. The 
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baryogenesis temperature is 



T B ~ Mi ~ 10 10 GeV , (17) 

and the washout factor k depends on the neutrino masses in a way, which will be 
discussed in detail in the follow ing chapter. Its determination requires the solution of 
the Boltzmann equationJ^^; in the estimate (JT5j) we have assumed a typical value, 
K ~ 0.1. The correct value of the baryon asymmetry, tjb ~ 10 -9 is then obtained as 
consequence of a large hierarchy of the heavy neutrino masses, which leads to a small 
CP asymmetry, and the kinematical factors / and The baryogenesis tempera- 

ture Tg ~ 10 10 GeV corresponds to the time ts ~ 10 -26 s, which characterizes the 
next relevant epoch before recombination, nucleosynthesis and electroweak transition. 

An important question concerns the relation between leptog enesis and neu trino 
mass matrices which can account for low-energy neutrino datcPElESEHEH . At 
present we know two mass differences for the light neutrinos and we have some in- 
formation about elements of the mixing matrix U in the leptonic charged current. 
Since U could be entirely due to mixings among the charged leptons, this does not 
constrain the light neutrino mass matrix in a model independent way. The light neu- 
trino masses can be either quasi-degenerate or hierarchical, and they can easily be 
consistent with successfull leptogenesis. In grand unified theories, due to quark-lepton 
unification, the hierarchical quark masses together with the seesaw relation require 
also hierarchical heavy Majorana masses. Neglecting the mixing of charged leptons, 
the constraints become very stringent and successful leptogenesis is possible only for 
a very small part of parameter spacd^ . In this case the enhancement of the CP 
asymmetry for partially degenerate heavy neutrino cP^^^ plays an important role. 
On the other hand, in unified theories large mixings between charged leptons easily 
occur. For instance, in a six- dimensional SO(10) model the parameters used in the 
above estimates for the baryon asymmetry were recently obtained^!. 

3. Quantitative Analysis and Bounds on Neutrino Masses 

Leptogenesis is a non-equilibrium process which takes place at temperatures 
T ~ Mi. For a decay width small compared to the Hubble parameter, T\(T) < H(T), 
heavy neutrinos are out of thermal equilibrium, otherwise they are in thermal equi- 
librium. A rough estimate of the borderline between the two regimes is given by 
Ti = H(Mi)3*. This is equivalent to the condition that the effective neutrino mass, 

m = (m ^ )n , (18) 
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equals the 'equilibrium neutrino mass' 



16vr 5 / 2 , ,„ tr 



.2 

m* = ^9 l J 2 TT- ~ 10_3 eV • ( 19 ) 



3v/5 M pi 

Here we have used M p i = 1.2 x 10 19 GeV and = 434/4 as effective number of 
degrees of freedom. For rh\ > m* ( rh\ < m*) the heavy neutrinos of type Ni are in 
(out of) thermal equilibrium at T = Mi. 

It is very remarkable that the equilibrium neutrino mass m* is close to the neutrino 
masses suggested by neutrino oscillations, a/ Am^ ol ~ 8 x 10~ 3 eV and ■sj Am 2 tm ~ 
5 x 10~ 2 eV. This suggests that it may be possible to understand the cosmological 
baryon asymmetry via leptogenesis as a process close to thermal equilibrium. Ideally, 
AL = 1 and AL = 2 processes would be strong enough at temperatures above Mi 
to keep the heavy neutrinos in thermal equilibrium and weak enough to allow the 
generation of an asymmetry at temperatures below M\. 

In general, the generated baryon asymmetry is the result of a competition be- 
tween production processes and washout processes which tend to erase any gen- 
erated asymmetry. Unless the heavy Majorana neutrinos are partially degenerate, 
M 2i 3 — Mi <C Mi, the dominant processes are decays and inverse decays of Ni and 
the usual off-shell AL = 1 and AL = 2 scatterings. The Boltzmann equations for 
leptogenesis read, 



dN 



dz 
<1Nb-l 
dz 



-(D + S)(N Nl -N%), (20) 
-si D (N Nl - N C N \) - W Nb-l • (21) 



Here Ni are number densities, z = Mi/T and D/(Hz), S/(Hz) and W/(Hz) denote 
decay rate, AL = 1 scattering rate and AL = 2 washout rate, respectively; all rates 
are normalized to the Hubble parameteJ^ . 

In order to understand the dependence of the solutions on the neutrino parameters, 
it is crucial to note that the rates D/(Hz), S/(Hz) and the resonance contribution 
(W — AW) /(Hz), are all proportional to fh\. One findiPJ, 

D,S,W-W«^, AW c M " M r 2 . (22) 
Here fhi is the effective neutrino mass (|18jl. and m is a quadratic mean, 



m 



2 — tr (m^rriv) = m\ + m 2 + m 2 . (23) 



Eq. (J22|) implies that, as long as AW can be neglected, the generated lepton asym- 
metry is independent of Mi. For quasi-degenerate neutrinos, with increasing m, the 
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figure 1a 
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Figure 3: Reaction rates in the thermal plasma compared with the Hubble parameter H as function 
of z = T '/Mi. The neutrino parameters are M\ = 10 10 GeV, fh\ — 10~ 3 eV. The two branches for 
Tw at small z correspond to the upper (lower) bounds T^y (T^). 

washout rate AW becomes important and eventually prevents successful leptogenesis. 
This leads to the upper bound on the absolute neutrino mass scale discussed below. 

The decay, scattering and washout rates are shown in fig. El as functions of z for 
a typical set of neutrino parameters, Mi = 10 10 GeV, m% = 10~ 3 eV, fa = 0.05 eV. 
All rates are of order the Hubble parameter at z ~ 1 where baryogenesis takes place. 
The generation of the B — L asymmetry for these parameters is shown in fig. 0] for 
|ei| = 10~ 6 and for two different initial conditions: zero and thermal N\ abundance. 
The figure demonstrates that the Yukawa interactions are strong enough to bring 
the heavy neutrinos into thermal equilibrium before leptogenesis takes place. The 
resulting asymmetry is in accord with observation, t]b ~ 0.01 x Nb-l ~ 10~ 9 . 

Gi ven th e heavy neutrino mass Mi, the CP asymmetry e\ satisfies an upper 
bonne™. Using this bound one can determine the maximal baryon asymmetry 
rfg** as function of the masses fhi, Mi, and rrP^, 

Vb < ^ ax (mi,Mi,m) ~ 0.96 x 10~ 2 e^fhu M u m)K(fh u Mi, m) . (24) 

Requiring the maximal baryon asymmetry to be larger than the observed one, 

v ^(mi,Mi,m)>^ MB , (25) 
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figure 1b 




z=M7T 



Figure 4: Evolution of the N\ abundance and the B — L asymmetry for e\ = —10 6 and m = 0.05 eV. 
The full (dashed) line corresponds to zero (thermal) initial N± abundance. 



yields a constraint on the neutrino mass parameters mi, M\ and m. 

The maximal CP asymmetry as function of rhi, Mi and m is given by^^ 



3 Mim 3 

167T V 2 



. mi 
m 3 



1 + 



2 2\ 1/2 

m 3 — m l 

~ 2 

mi 



(26) 



For neutrino masses with normal hierarchy one has 



ITLn 



Am 



atm ' 



m 



777.1 



(27) 



and the dependence on m reads 



mo 



m 9 



- (m 2 + 2Aml tm + Am^,) , 



(m 2 - Aml tm + Am 2 sol ) , 
(™ 2 -A™L-2Am s 2 ol ) • 



(28) 
(29) 
(30) 



Note that E\ = for mi = mi. In general, one expects m± < mi < m%. Here the 
lower bound ' ' is always true whereas the upper bounds holds if there are no strong 
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Figure 5: Normal hierarchy case. Curves, in the (mi-Mi)-plane, of constant r/g ax = 10 -10 
(thin lines) and r/^ ax = 3.6 x 10~ 10 (thick lines) for the indicated values of fn. The filled 
regions for ry^ ax > 3.6 x 10~ 10 are the allowed regions from CMB. There is no allowed region 
for m = 0.20 eV. 

cancellations due to phase relations between different elements of the neutrino mass 
matrix. 

Using the upper bound ()26|) on the CP asymmetry one can calculate the maximal 
baryon asymmetry. The CMB constraint (|25|) yields for each value of ma domain 
in the (m 1 -M 1 )-plane which is allowed by successful baryogenesis (cf. fig. EJ). For 
m < 0.20 eV this domain shrinks to zero. Using the relations (j28|) - (p?U|) . one can 
easily translate this bound into upper limits on the individual neutrino masses, 

mx,m 2 < 0.11 eV , m 3 <0.12eV. (31) 

Note that these bounds are a factor of factor of two below the recent upper bound 
of 0.23 eV obtained by WMAF . For an inverted hierarchy of neutrino masses one 
finds very similar upper bounds, m\ < 0.11 eV and 1712, < 0.12 eV. In a complete 
analysis the change of neutrino masses between the mass scale of leptogenesis and 
the electroweak scale has to be included. Generically, one expects that this will make 
the upper bounds on the neutrino masses more stringentPSJ . 

In a similar way one can obtain a lower bound on M\, the smallest mass of the 
heavy Majorana neutrinos. One findJ^ 

Mi > 4 x 10 8 GeV . (32) 
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Figure 6: Washout factors as function of Zi — Mi/Tj without (full line) and with (dashed line) 
iVi-top scatterings; Mi = 10 10 GeV. 

As a consequence, thermal leptogenesis requires a rather high reheating temperature, 
T R > T B ~ Ml 

An important question for leptogenesis, and baryogenesis in general, is the depen- 
dence on initial conditions. This includes the dependence on the initial abundance 
of heavy Majorana neutrinos, as discussed above, and also the effect of an initial 
asymmetry which may have been generated by some other mechanism. It turns out 
that heavy Majorana neutrinos can efficiently erase initial asymmetries. 

For large initial asymmetries, one can neglect the small asymmetry generated 
through the CP violating interactions of the heavy neutrinos, i.e. one may set e± = 0. 
The kinetic equation (|2*Tj) for the asymmetry then becomes 

^ = -WN B . Lt (33) 

where —Nb-l is the number of lepton doublets per comoving volume. The final B — L 
asymmetry is then given by 

N B _ L = 0( Zi )N B _ L , (34) 

with the washout factor 

0(z i ) = e-^ dzW(z) . (35) 
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The result of a quantitative analysis^* is shown in fig. The washout becomes 
very efficient for rh\ > m* ~ 10~ 3 eV. Already for fhi = 5 x eV one has 

w(zi) ~ 1CT 7 . . . 1CT 5 , indicated by the dashed and full lines, respectively Hence, an 
initial asymmetry several orders of magnitude larger than the presently observed one 
can be reduced to a value below the one generated in leptogenesis. The range for 
w(z{) is due to a theoretical uncertainty in the treatment of iVi-top scatterings. Note, 
that a plateau for w{z\) is reached for values of Z\ just below one. 

We conclude that leptogenesis naturally explains the observed baryon asymmetry 
for neutrino masses in the range 

1(T 3 eV < rrn < 0.1 eV , (36) 

almost independent of possible other pre-existing asymmetries. It is very remarkable 
that the data on solar and atmospheric neutrinos indicate neutrino masses precisely 
in this range. 
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